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Abstract 

A theory has been developed to calculate the Hadamard- Walsh transform from a list 
of cubes specification of incompletely specified Boolean functions An efficient algorithm 
to generate disjoint cubes from non-disjoint ones of a cover has been developed The 
transformation algorithm makes use of properties of list of disjoint cubes and allows the 
determination of the spectral coefficients in an independent way 

Decomposition of Boolean function at the functional level is a difficult problem since 
it requires a global approach We discuss linear decompositions of Boolean functions imple- 
mented as Programmable Logic Arrays (PLAs) Results with experimentation with certain 
functions shows that functions with imbedded addition are most benefited from this ap- 
proach 

Syndrome testing is particularly useful in two level circuits such as PLAs It is 
shown that weighted sum of syndromes of all the outputs covers all single stuck-at-faults, 
bridging faults and cross-point faults An algorithm is presented, which checks for the kps'.? 
testability of a PLA and suggests hadware modifications, if found WS'5-untestable 
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Chapter 1 


INTRODUCTION 


During the past three decades a new field of digital logic theory has been born and its 
boundaries continuously extended There was an eager attempt to apply orthogonal trans- 
form techniques to the design, optimization, and testing of digital logic It was conjectured 
that these so called spectral methods would provide a unified approach to the synthesis of 
anlog and digital (binary and multivalued) circuits, making it possible to apply the power- 
ful Fourier transform techmques, widely used m analog circuit synthesis, to digital network 
design Two motives for this work may be discerned first, the purely academic one of 
applying existing mathematical techmques to the field of digital logic design, and second, 
the increasing appreciation of the limitations of existing algebraic and geometric methods 
m handling digital data for logic network design purposes The theory that subsequently 
emerged has proved to be fascinating 

Digital logic design is a peculiar discipline in that the logic design for a given re- 
quirement, which may be in total a very large system, is frequently made with the use of 
sophisticated design procedure Very great sophistication, however, may be present in the 
logical verification of the complete assembly and m its translation into a microelectronic 
realisation Though the increasing capabilities of LSI and VLSI fabrication have so far kept 
pace with the designers’ need, there has always been a rising demand for ’’better” design 
techmques 

However, the spectral methods of digital logic theory represents a modern approach 
to expressing conventional digital data, one which can provide various insights into the 
structure of the data which are absent from the classical Boolean algebra and truth-table 
formats The pioneering work, particularly that of Karpovsky and Lechner, forms the 
basis of this approach Others have contributed and amplified the basic concepts and have 
translated the under lying mathematics into engineering tools which may be more acceptable 
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to a digital logic designer 

Spectral techniques have been applied to Boolean function classification [1], [2], [3], 

[4] , disjoint decomposition [2], [5], parallel and serial Lnear decomposition [2], [6], [7], [8], 

[5] , multiplexer synthesis [2], [9], prime imphcant extraction by spectral summation [2], 
[8], [10], threshold logic synthesis [1], [8] and state assignment [7], [11] Spectral methods 
for testing of logical networks by verification of the coefficients in the spectrum have been 
developed [2], [6], [8], [12], [13], [14], [15], [16], [17], [18], [l9] It should be stressed that 
an important problem of finding the complement of a Boolean function that has high com- 
plexity in the Boolean domain [20], [21] can be solved very easily in the spectral domain 
because complementing the Boolean function corresponds to changing the sign of every 
spectral coefficient Tautology of a Boolean function can be verified by calculating a certain 
coefficient (DC coefficient) The renewed interest in applications of spectral methods in 
logic synthesis is caused by their excellent design for testability properties and the possibil- 
ity of performing the decomposition with gates other than the ones used in most classical 
approaches Spectral techniques have also been used for data transmission, digital filtering, 
image processing and pattern analysis This thesis will attempt to introduce the underlying 
theory of this area, that of orthogonal transforms and resulting spectral data While the 
underlying theory of this area is apphcable to any-valued logic, and not exclusively to the 
binary case, we will confine our discussion herewith to two- valued digital netwotks 

The second chapter gives a complete explanation of the spectra of discrete functions 
and their calculations An efficient algorithm for calculation of Hadamard-Walsh trans- 
form from a cube array specification of incompletely specified Boolean functions has been 
designed 

In the third chapter we discuss the linear decomposition of Boolean functions at 
the functional level with emphasis to Programmable Logic Arrays (PLAs) An procedure 
for the hnear decompositon is presented, of Boolean functions Experimental results are 
presented for some commonly used functions implemented as PLAs which estabhsh that 
such decomposition can often result in imroved implementation of logic functions 

The most recent appbcation for spectral techmques has been in fault diagnosis This 
exciting area is dealt with in the fourth chapter One of the serious problems with digital 
circuits concerns their testing There are two aspects to this - first to verify whether or 
not a circuit is performing correctly and, second, if there is a fault, to find its location 
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(fault isolation). Of course, from the user’s point of view there is not usually much purpose 
in locating a fault beyond identifying a chip that needs replacing There is no need to 
identify the the exact location of the fault in the chip We shall discuss about a simple and 
effective fault detection technique known in the bterature as Syndrome testing We apply 
this technique with the algorithms discussed in second chapter to fault detection of PLAs 
All the algorithms discussed in the second, third and fourth chapters are suitable for 
automation and have been implemented in ‘C’ language 
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Chapter 2 


THE SPECTRA OF DISCRETE FUNCTIONS 


In this chapter we shall review the theoretical aspects of spectra of discrete functions in 
general and later on formulate an algorithm to calculate the spectra of incompletely spec- 
ified Boolean functions represented either as a bst of true/false minterms or as an already 
minimized sum-of-products Boolean expression (SOPE) 

2.1 Boolean and Spectral domains 

The majority of existing methods for design and analysis of switching circuits are concerned 
with the properties of Boolean functions since it was proved that the Boolean domain 
provides a precise model for the analysis of switching circuits The behaviour of a device is 
represented by a function !F{xi^X2^ ,^n) of its input variables a:i,a:2, This function 

can most convemently be defined by a truth table For example, a function ^(xi, X2^ 2^3, X4) 
of four variables xi^X2^xs and X4 is illustrated in Fig 2 1 

It is common to use the product and sum operators of the Boolean algebra together 
with negation to define such functions-for example, J^{xi^X2) = X1X2 + ^12:2 Part of 
the problem with the defimtion in the Boolean domain is that each of the entries in the 
column for /■(a:i,X2,a;3,a:4) in Fig 2 1 teUs us precisely the behaviour of the function 
at a single point but nothing of its behaviour at other points It is possible to give an 
alternate representation of a function where the information about the function is much 
more global in nature This alternate representation is in the spectral domain, and it 
has been demonstrated [2], [7] that a number of properties are much easily deduced in the 
spectral domain than in the Boolean one 

The basic idea of the spectral domain, and how to get there, is illustrated in Fig 
2 2 K we are to avoid losing information, we shall have to ensure that the transform can 
be reversed, 1 e , that we can move to and from the spectral domain without any loss of 
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Figure 2 1 The function 0:2, 2:3, 0:4) 


information 

The question arises as to the reason for considering the spectral domain and if there 
are any real purposes for its use To understand the first difference between the Boolean 
and spectral domain, let us consider a Boolean function T{X) of n variables One row 
of the table defining this function provides complete and precise information about the 
behaviour of the function for one combination of the input variables Of course, it does 
not tell us anything about the value of the function anywhere else The combination of 
the knowledge of the behaviour of for the 2 ^ rows of the table gives a complete 

definition of the functions Similarly the spectrum for a function of n variables also contain 
2” values, which together completely define the function and can be used to recover its 
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Figure 2 2 The transform 
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Boolean specification Each of the 2^ values in the spectrum (the spectraj coefficients) 
contains some information about the behaviour of the function at all 2^ points, but does 
not contain complete information about any of them The combination of all the values 
in the spectrum does lead to complete information about the function, but each individual 
coefficient gives us some global information about the whole function 

2.2 The Transformation of Binary Data 

As shown in Fig 2 1 any given combinatorial function JF(X) may be explicitly defined by 
its truth table, which lists all 2^ combinations, against each combination being given the 
local function output value X(X) Conventionally, all a:, input variables and the output 
value of J^(X) are expressed in {0,1} notation Note that each input combination and the 
corresponding output value of X(X) is discrete data, no information about the output value 
of at any input combination being obtaanable from any other input combination 

With the order of tabulation of the input variables standardised in straight binary 
code (SBC), the 2^ local output values may be treated as a column vector Z defimng 
X(X) such that the conventional { 0, 1, — } values (false, true and don’t care mmterms) 
correspond to { 0, 1,0 5 } codings respectively as in Fig 2 3(a) Should the binary logic 
values be recoded from { 0, 1, — } to { +1,— 1,0 }, respectively, then an alternative column 
vector Y defimng ^(X) is obtained, as shown in Fig 2 3(b) The reason for such recoding 
can be found in [2], [7] 

In both Y andZPwe have 2’^ntries necessary to fully define an n- van able function 
J^(X). As an altrnative to these 2^ efntries, we may also specify 2^ coefficients that define 
^(X), these coefficients being the spectral coefficients of X{X)^ these coefficients may also 
be bsted in a defined order to give a column vector representation, which is identified by R 
or S As will be seen, the entries in R and S will no longer be binary values as in the local 
value truth table vectors Y and Z 

The defimtion of the individual spectral coefficients in R and S for any J^(X) may 
be approached from several mathematical viewpoints Here we shall first consider the 
Hadamard transform T’^, which will directly generate S from Y or R from Z, following 
which we consider alternative but equivalent mathematical definitions and aspects 
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Figure 2 3 Truthtable formats for example four-variable function 
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2.2.1 The Hadamard Orthogonal Transform Matrix 

The Hadamard trajisform is a complete orthogonal square matrix, with row and column 
entries € {+1,-1}, and with a recursive structure as follows 


T 


n __ 


jjin— 1 _'pn— 1 


(2 1 ) 


Note that the dimensions of the transform are 2^X2^^ for any n For increasing n we have 
the real integer values 


T° = +1 



1 - 1-1 1 


We may alternatively express this recursive structure by 



(2 2 ) 


(2 3) 


where (gi here denotes the Kronecker product operator 

The value of each element tjk m the Hadamard matrix may be individually defined 
by mathematically summing the element-by-element multiplication of the binary (mod 2) 
expansions of j and k (call this v) and taking the uth power of —1 For example, taking 
the third row (y = 10), fourth column (k = 11) entry <23 of T", shown below we have 


k{binary) 00 01 10 11 


j (binary) 


00 

■ 1 1 1 1 ■ 

01 

1-1 1-1 

10 

1 1 - 1-1 

11 

1 - 1-1 1 

11 , 1 ) = { 1 1 - 1 - 0 1 } = 
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This may formally be expressed by 


tjk = 

(2 4) 

where determined by the binary expansions of j^k respectively, y, A: 6 0 to 2’^«1, 

and where 


+ Jo2°} 

k = {^n-l2" ^ fc„_22" ^ -f 

+ fco2°} 


The individual row vectors of T” define a complete orthogonal series of functions, 
which IS closed under row multiplication The inner product of any two matrix rows jajp 
has the orthogonal property 


2^-1 

Ar=0 


2^ if a = /3 
0 otherwise 


(2 5) 


This inner product relationship holds also for any two- column vectors 

The Hadamard transform matrix has the additional property that its transpose is 
identical to itself, that is, T^ = [T’^]^ thus giving the simple inverse transform relationship 
for any n of 

IT")-' = ^[T”) (2 6) 

2.2.2 The Hadamard Transformation of Binary Data 

The Hadamard matrix may be used to transform any 2" column vector of numbers into an 
alternative data column vector No information is lost m this transformation, the original- 
data being fully recoverable from the resulting vector by application of the inverse transform 
[T"]~i We will use it in the transformation of a column vector Y or Z representing the 
output truthtable of a given function 

Example Considering the four-variable example function of Fig 2 3 and taking the 
(conventional) 0, 1 coding we may transform Z with an n = 4 Hadamard transform as 
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follows 
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Now, 


R‘ = [T‘‘zf 

= [l2 5 -1 5 -3 5 -1 5 -15 0 5 -1 5 0 5 

05 -1 5 05 -1 5 25 05 25 05] 


(2 7) 

(2 8 ) 


If tlie output vector of is recoded from { 0, 1, — } to { +1,-1, 0 } then we have, 


Y* 


0 -1 -1 -1 1 - 1 - 1-111 -1 -1 -1 -1 -1 -1 ] 


S‘ = [T'‘Y]‘ (2 9) 

= -9 3 7 3 3 -1 3 -1 -1 3 -1 3 -5 -1 -5 -1 ] (2 10) 

The resulting vectors spectra R and S each uniquely represent the given function T{X) 
The individual coefficient values in R and S are linearly related due to the linear recoding 
between Z and Y, the relationships being 


Vj = 

l-2zj 

(2 11) 

To = 

|(2" - »o) 

(2 12) 

Toe = 

1 

(2 13) 
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The inverse transformation back from the spectraldomain to the Boolean domain is directly 
given by, 


Z = [T"]~^R 


(2 14) 


[Tnj-lg 




-It 


(2 15) 


Whilst the coefRaents R or S of equations ( 2 8) or ( 2 10) uniquely represent a single 
combinatorial function ^(X), a joint spectrum of two or more functions may be computed 
Example Let us consider the two functions J^i(X) and E 2 {X) tabulated below En- 
coding and combimng the two individual outputs we have, 
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The Hadamard transform now gives 
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(2 16) 


ij»n 


n+2 


R 


1+2 


2.2.3 The Fast Transform Procedures 


The step-by-step execution of Hadamard forward or inverse transform involves the ± sum- 
mation of a total of 2^X2”' individual product terms Thus the computation of spectra using 
matrix operations is impractical for all but extremely small functions Due to the recursive 
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structure of [T], however, it is possible to use a fst,^^ place transformation algorithm [7], 
[ 22 ], as given below 


to{u>) = lF{u>) 

(2 17) 

to( 2 "-^-l-a;) = E{2^-^+u) 

(2 18) 

= t,_i( 2 a;) -b tt-i( 2 u) -f 1 ) 

(2 19) 

t,( 2 " ^- 1 - 0 )) = t,_i( 2 w) — t,_i( 2 w -h 1 ) 

(2 20 ) 

a; = 0 , 1 , , 2 ""^-l 

(2 21 ) 

1 = 1 , ,n 

(2 22 ) 


On completion the array t contains the transform Instead of involving the total number of 
individual product terms, the total number of operations is reduced from 2"X2” to 2’^Xn, 
which IS expressed as N/ 052 N, where N = 2 " 

Example This is illustrated below for the function E{X) of Fig 2 3 for S spectrum 



H^) 

to{u;) 
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-1 

-3 

-5 

-9 
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1 

-1 

-1 

-2 

-2 

-4 

3 

Si 
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-1 

-1 


0 

3 

7 

S2 
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-1 

-1 

-2 

-4 

0 

3 

S 12 

4 

1 

1 

2 

1 

3 

3 

S 3 

5 

-1 

-1 

-2 

2 

4 

-1 

Sl3 
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-1 

-1 

-2 

0 

3 

3 

S 23 

7 

-1 

-1 

-2 

0 

0 

-1 

Sl23 
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1 

1 

1 

1 

-1 

-1 

S 4 

9 

1 

1 

0 

2 

4 

3 

Sl4 

10 

-1 

-1 

2 

4 

-1 

-1 

S 24 

11 

-1 

-1 

0 

0 

0 

3 

Sl24 

12 

-1 

-1 

0 

1 

-1 

-5 

S 34 

13 

-1 

-1 

0 

2 

4 

-1 

Sl34 

14 

-1 

-1 

0 

0 

-1 

-5 

S 234 

15 

-1 

-1 

0 

0 

0 

-1 

<*1234 


Thus the obtained spectrum S = ^ 4 ( 0 ?) is exactly the same as the one obtained in the 
previous case 

2.2.4 Properties of S and R Spectral Coefficients 

The principal properties of the coefficients of the S and R spectrum for completely and 
incompletely specified Boolean functions for the well-known Walsh-type transform matrices 
are shown below according to [1], [2], [7], [ 8 ], [23], [25] 
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1 The transformation matrix for each ordering of the Walsh functions is complete and 
orthogonal^ therefore, there is no information lost in the S and R spectra concerning 
the minterms of the Boolean function T 

2 Only the Hadamard- Walsh matrix describing the Hadamard-Walsh transform hcis the 
recursive Kronecker product structure Other possible variants of the Walsh trans- 
forms, described by the corresponding matrices, are known in the literature as the 
Walsh-Kaczmarz, Rademacher-Walsh, and Walsh-PaJey transforms 

3 Only the Rademacher-Walsh matrix is not symmetric^ all other variants of Walsh 
matrices are symmetric, so that, disregarding a scaling factor, the same matrix can 
be used for both the forward and inverse transform operations 

4 Each spectral coefficient sj or r/ is described by its order^ subindexes and magnitude 
The order of the spectral coefficient is equal to the number of subindices, and the 
subindices are the subscripts of all variables of a standard trivial function correspond- 
ing to the coefficient The magmtude of a spectral coefficient is its value. In this 
chapter, i, j and k denote subindices, and the order is denoted by o 

5 When the classical matrix multiplication method is used to generate the spectral 
coefficients for different Walsh transforms (different T matrices represejit^ different 
Walsh functions with different orderings), the only difference in “^he rtesult is th4 
ordering of the spectral coefficients The coded vector Z or Y corresponding to the 
original truth vector of a Boolean function is the same for all orderings of Walsh 
functions The values of the coefficients sj and r/ with the same subindices are the 
same for every ordering of Walsh transforms 

6 Each spectral coefficient (either sj or r/) gives a correlation value between the Boolean 
function T and a standard trivial function uj corresponding to the coefficient The 
standard trivial functions for the spectral coefficients are respectively 

• for the dc coefficients (direct current coefficients) sq or ro, the universe of the 
boolean function denoted by uq 

• for the first order coefficients sj or rj{I = 2 ,z ^ 0), the variable of the boolean 
function !F denoted by u, 
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• for the second order coefficients sj or rj{I = tj, t 0, j ^ 0, t / j), the Exclusive- 
OR function between variables x, and Xj of the boolean function ^ denoted by 

• for the third order coefficients sj or r/(/ = tjk^t ^ 0,jf 0, A: 0,z i 7^ 

^ fc), the Exclusive-OR function between variables x,, Xj and Xk of the 
boolean function denoted by ujA, etc 


7 The number of spectral coeficients of 2:th order is equal to 
the number of variables of a Boolean function 


n 

z 


where n is 


8 For a completely specified Boolean function the maximal value of any individual spec- 
tral coefficient sj is +2”^ and the minimal value is —2'^ These values occur when the 
Boolean function T is equal to either a standard trivial function uj (sign +) for the 
maximal value or to its complement (sign -) for the minimal one In either case, all 
the remaining spectral coefficients have zero values because of the orthogonality of 
the transform matrix T 


9 The maximal value of any spectral coefficient rj except tq is +2’^^^ and will residt 
when the Boolean function !F is equal to complement of a standard trivial function 
uj The mimmal value is and will result when T is equal to a standard trivial 

function uj In either case, all the remaining spectral coefficients have zero values 
because of orthogonality of the transform matrix T 

10 For an incompletely specified Boolean function the maximal value of any individual 
spectral coefficient sj is +2” — 1, and the mimmal one is —2’^ + 1 These occur when 
the Boolean function ^ is equal to either a standard trivial function for the maximal 
value of sj or to its complement for the minimal value, in all but one iUlnterm which 
is a don’t care 


11. The maximal value of tq spectral coefficient is +2"^ and will result when the Boolean 
function .7^ is a tautology The minimal value is —2’^ and will result when !F is equal 
to the complement of the tautology The tautology is the logical function for which all 
the imnterms are true 


12 When, for more than half of the spectraol coefficients of any completely specified 
Boolean function .7^, the majority of the mmterms have the same logical values as the 
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minterms of standard trivial functions, the sum of all the coefficients of the S spectrum 
has a maximum value equal to +2” When, for more than half of the spectraol 
coefficients, the majority of the minterms of T have the complemented logical values 
to the minterms of standard trivial functions, the sum of all the coefficients of the S 
spectrum has a mimmal value equal to 

13. For any incompletely specified Boolean function, the sum of all coefficients of the S 
spectrum has a maximal value of +2” — 1 and a minimal value of —2’^ + 1 The 
maximal or mimmal value happens when the Boolean function has exactly one don’t 
care minterm and all the spectral coefficients follow the rule of property 12 

14. With the exception of uq, every standard trivial function uj corresponding to an n 
variable Boolean function !F h 2 LS the same number of true and false minterms That 
number is equal to 

15. For each true minterm the coefficients from the spectrum S are Sq = 2’^ — 2, and all 
remaimng 2^^ ~ 1 spectral coefficients s/ are equal to ± 2 The method for choosing 
the signs of spectral coefficients is described in Section 2 4 

16. The spectrum S of each false minterm is given by s/ = 0 

17. For each don’t care minterm the coefficients from the spectrum S are sq = 2’^”^ — 1, 
and all remaining 2^^ — 1 spectral coefficients 6/ are equal to ± 1 The method for 
choosing the signs of spectral coefficients is described in Section 2 4 

The new method of calculating spectral coefficients of Boolean functions described 
in the Section 2 4 needs the representation of Boolean functions in the form of arrays of 
disjoint cubes The method gives the correct values of spectral coefficients independently 
from the shape and size of disjoint cubes in the array of cubes as long as all the minterms 
of the function are covered only once by a cube The next section describes the algorithm 
that generates this kind of representation of Boolean functions originally represented either 
as a list of true/false minterms or as an already minimized SOPE But before we describe 
the algorithm, here are some notations and defimtions [26] 

Definition 1 A cube is a set C of bterals such that x £ C impbes that x is not an 
element of C {e.g {a,6,c} is a cube but {u, a} is not a cube) A cube represents the 
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conjuction of its literals. In a cube^ bterals 0 and 1 are referred to as bound coordinates and 
X (don’t care) as free 

Definition 2 A ON-cube is a cube such that the Boolean function which consists of 
this cube evaluates to 1 for the combination of the input variables represented by this cube 
Similarly we have OFF-cubes and DC~cubes 

Definition 3 A vertex is a special cube in which all coordinates (literals) are 0 or 1 
Definition 4 A cube is Ca is said to contain another cube Cfe, if Cb can be obtained 
from Ca by appropriately changing all don’t cares (represented here by x) in the cube Ca 
to Is or Os (e g lOx contains both the cubes 101 and 100) Note that, both the cubes Ca 
and Cb should be of the same type {i e ON, OFF or DC) 

Definition 5 The intersection of cubes Ca and Cb is defined as the smallest cube 
containing all of their common vertices - it is contained in all cubes containing them If they 
have no common vertices, they are said to be disjoint Using the symbol I for intersection, 
{lxl)I(0xl) = <!> Thus the cubes 1x1 and Olx are disjoint Similarly, (a;10a:)/(01a:a:) = 
010a: 

Definition 6 The sharp product (termed as ^-product) of cube Ca with Ch is a 
differencing product, consisting of a set of vertices of Ca which are not contained in Cb 
Again we assume that both the cubes are of the same type There are three rules 

• If Ca and Cb are disjoint, then Ca#Cb = Ca, 

• If Ca IS contained in Cfc, then Ca#Cb contains no cubes, i e Ca#Cb = 0 

• If neither of the above conditions hold, then the output will be a set of certain faces of 
Ca that are not in Cb For each coordinate i for which Ca{t) = x and Cfe(t) = 1 or 0, 
there will be a cube Cd m the ^-product having Cd{i) — e and all other coordinates 
identical to those of Ca 

Example 


(lxlx)#(xlx0) = lOlx, 1x11, 

(xxx):jl^(101) = Oxx, xlx, xxO 

Definition 7 The disjoint sharp operation (termed as -product) aids m the counting 
the number of vertices When is formed between two cubes Ca and Cb (of same type) it 
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obtains a set of disjoint cubes. To effect Ca#jC6, first the ordinary Ca#Cb is formed Let 
the first cube Ci produced be as with the ^^t-product, the second cube C2 is modified so that 
in the first coordinate t(l) where C\ has a 1 (or 0) and C2 has an x, this x is changed to 
an opposite 0 (or 1) in the product For the third cube C3 produced by the ^-product, 
the t(l) coordinate is made similarly equal to 0 (or 1), likewise the second coordinate t{2) 
IS found where C2 has a 1 (or 0) and C3 has an x This x is changed to a 0 (or 1), rendering 
it disjoint from Ci and C2 
Example 


xxxx 

xxxx 

# 1011 #, 

1011 

Oxxx 

Qxxx 

xlxx 

liar a: 

xxOx 

100a: 

xxxO 

1010 


2.3 Algorithm to Generate Disjoint Cubes 

As Will be seen the calculation time for Hadamard spectrum increases linearly with the 
number of disjoint cubes Thus, the determination of a quasi-minimal number of cubes in 
the disjoint cube representation of a Boolean function is crucial for the effective calculation 
of its spectrum 

Several algorithms for the generation of arrays of disjoint ON- and DC-cubes (if 
any) or an array of disjoint OFF-cubes were described and used in ESPRESSO [20], [24], 
EXORCISM [27] and PALMINI[28] Here the algorithm and its implementation as described 
in [24] are improved The drawbacks of this algorithm are also overcome. 

A logical function is represented in the arrays of cubes form It is enough to have only 
any two of the three possible sets of cubes (ON, OFF and DC) In the following presentation 
of the algorithm it will be assumed that only the ON- and DC-arrays are available (the 
same as assumed m [20]) This algorithm can be applied to the ON-array and DC-array 
separately maintained as umdirectional hnked hsts 

The random ordering of cubes is replaced by a list of cubes in which the are sorted 
according to their size Every time two cubes Ca a.nd Ch are compared The algorithm uses 
two pointers a and h which correspond to the actual position of cubes Ca and Ch m the hst 
The pointer a indicates the position of the cube which is the first of the pair of the cubes 
being compared, the pointer b indicates the position of the second cube m the pair The 
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pointer a changes its values till the (n - l)th node (n = no of nodes in the list) is reached 
and b changes its values from 2nd to nth node accordingly At a given moment, the value 
of the pointer b is always greater by at least one than the value of the pointer a During the 
execution of the algorithm the number of nodes in the list can be different from the original 
number of cubes due to the following reasons 

• the disjoint sharp operation -operation) can generate more than one cube as its 
result 


• the contain operation can remove some cubes 


As a last step of the algorithm, the cubes are merged, where possible, to obtain a smaller 
total number of disjoint cubes 

Algorithm Generation of a set of disjoint cubes from the set of non-disjoint cubes 


initialise 

a = head of the list, 

if(a == end-of-bst) stop, /* if the list is empty 
b = a next, 

if(6 == end-of-list) stop, /* if the list has only 
do { 

d = step = 1, 
count = 0, 
do { 

if Ca and Ch are not disjoint 

{ 

if Ca contains Ch 

{ 

remove Ch from the list and 
update the list, 

} 

else 


V 

one node */ 


{ 

do C^'^jCh (2 

d = no of cubes resulted from Ctp#jCh>^^ 

i{{count == 0) step = d, 

replace Ch by the first solution cube. 

Insert the other solution cubes in the list to the next of C^, 
6 = 1,/* advance b so that it points to 

the node containing the last solution cube */ 


} 

} 

b = b next, 
count + +, 

} while (b NULL), /* until the last node is processed */ 
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b = a + step + 1 , 
a = a —y next, 

} while (a end-of-list), 
merge the cubes, 

Example The above algorithm is applied to the function of Fig 2 3^ The steps of 
the execution are shown in Fig 2 4 Fig 2 4(a) shows the minimized function of Fig 2 3 
minimized by ESPRESSO [20] A in the hst of cubes indicates that the disjoint sharp 
operation (#_,) has to be perfomed between those two cubes Fig 2 4(e) is the hst after 
the disjoint operation Fig 2 4(f) is the hst after merging the possible cubes 
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ON 

*00x1 

ON 

0000 

DC 
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xxlx 

ON 

IlOx 

ON 

0x01 

ON 

0000 
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Figure 2 4 Stages of execution of the algorithm to generate a disjoint cube representation 


2.4 New Method for Calculation of Spectrum 

An algorithm already exists for calculating spectral coefficients for completely specified 
Boolean functions directly from a SOPE [2], [29] When the imphcants are not mutually 
disjoint, this algorithm requires an additional correction to calculate the exact values of 
spectral coefficients for minterms of Boolean function E that are included more than once 
in some imphcants But in the disjoint cube representation, the exact values of spectral 
coefficients can be calculated immediately without having to perform correction operations 
Here the extension of the algorithm to incompletely specified Boolean function is proposed 
Definition 1 The cube of degree m is a cube that has m literals that can be either 
in affirmation or negation (i e , m is equal to the sum of the number of zeroes and ones in 
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the description of a cube) Let p denote the number of x’s in the cube and let n denote the 
number of variables of a Boolean function Then, n = m + p 

Example The cube 1x00 is of degree 3 since three of the literals describing this cube 
are either in affirmation (xi) or negation (xa and X4) The cube does not depend on literal 

Definition 2 The partial spectral coefficient of an ON- or DC-cube with degree m of 
a Boolean function is equal to the value of the spectral coefficient that corresponds to the 
contribution of this cube to the full n-space spectrum of J- The number of partial spectral 
coefficients npsc describing the T is equal to the number of ON- and DC-cubes describing 
this function 

Example Considering Tables 2 1 and 2 2 representing this method of calculation for 
the same function of Fig 2 4, each row shows the partial spectral coefficients of either On- 
or DC-cubes of a Boolean function This function has four partial spectra, which is equal 
to the number of disjoint ON- and DC-cubes describing this function (npsc = 4) 


X1X2X3X4 
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4 
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0 

0 

0 

4 

4 

-4 
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ON 

12 

4 

-4 

0 

-4 

4 

0 

4 

0 

-4 

0 
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7 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 
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-9 

3 

7 
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-1 

3 

-1 

3 

3 

-1 
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Table 2 1. Table showing DC, 1st and 2nd order coefficients 


xiX2a;3X4 


5123 

5124 

5134 

5234 

51234 

xlxx 

ON 

0 

0 

0 

0 

0 

xOll 

ON 

0 

0 

0 

-4 

0 

10x0 
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0 

0 

0 

0000 

DC 

-1 

-1 

-1 

-1 

-1 



-1 

3 

-1 

-5 

-1 


Table 2 2 Table showing 3rd order coefficients 

Suppose the hsts of disjoint ON- and DC-cubes that fully define ^ are given Then 
each cube of degree m can be treated as a minterm within its particular reduced m-space 
of T The spectrum of each true nunterm is given by 50 = 2^ - 2, and all remaining 
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2" — 1 coefficients are equal to ±2 (property 15) Similarly, the spectrum of each don’t care 
minterm is given by sdcq = 2 **“^ — 1 , and all the remaining 2 ’‘ - 1 coeffiaents are equal 
to ±/,l (property 17) The symbols sdci denote spectral coefficients corresponding to 
DC-cubes (when / = 0, the symbol sdco denotes a DC spectral coefficient corresponding 
to a DC-cube) 

Cubes of degree m have the following properties 

1 The contribution of an ON— cube of degree m to the full n-space spectrum of function 
^(where n is the number of variables in the function calf) is represented as follows 

So tn full n-space = 2 " - 2 x 2 ^ (2 23) 

and 

s/ in full n— space = sj tn m — space X 2^ (2 24) 

where 1^0 

2 The contribution of an DC— cube of degree m to the full n-space spectrum of function 
T IS related as follows 

sdco m full n-space = 2'*“^ -2x2’’ (2 25) 

and 

Sdci n-space = Sdci ^nm — space X 2’’ (2 26) 

where / 7 ^ 0 When the above formulas are applied to minterms (i e for m = n and 
p = 0) they reduce to properties 15 and 17. The contribution of a DC-cube of degree 
m IS equal to one half of the contribution of an ON-cube that has the same degree m 

Equations 2 24 and 2 26 determine the absolute values of those partial spectral 
coefficients sj that are not equal to zero for a given cube Properties 3-4 determine 
the signs of the partial spectral coefficients, and whether some of them are equal to 
zero 

Example Considering Tables 2 1 and 2 2 again, the value of partial spectral coeffi- 
cient So, corresponding to the On-cube lOarO (n = 4,p = 1) is equal to 2 “* — 2 x 2 ^ = 12 
according to Equation 2 23 The absolute values of those partial spectral coefficients 
sj that are not equal to zero are calculated according to Equation 2 24 and are equal 
to 2 X 2^ = 4 
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3 If m a given cube the x, variable of a Boolean function is denoted by the symbol 
then all of the partial spectral coefiicients sj whose indices I contain the subindex t 
are equal to zero 

4 For an ON- or DC-cube of degree m the number of nonzero partial spectral coefficients 
IS equal to 2^“^, except for p = n — 1 when there is only one nonzero partial spectral 
coefficient 

Example Consider Tables 2 1 and 2 2 again The variable X3 is denoted by symbol 
X in the cube 10x0 So all partial spectral coefficients with subindex 3 are equal to 0 
Therefore, S3 = S13 = S23 = ^34 = ^123 = ^134 = ^234 = ’§1234 = 0 For this cube by 
property 4, the number of non-zero partial spectral coefficients is equal to 2^”^ = 8 

The cube xlxx has three variables denoted by x xi, X3 and X4 Therefore, by 
property 3 and 4, only the partial spectral coefficient S2 is non-zero 

The following properties describe the signs of each partial spectral coefficient s/, where 

7/0 

5 If in a given cube the x* variable of a Boolean function is an affirmation, then the sign 
of the corresponding first-ordert coefficient is positive, otherwise for a variable that is 
in negation, the sign of the corresponding first order coefficient is negative 

6 The signs of all even-order coefficients are given by multiplying the signs of the related 
first-order coefficients by -1 

7 The signs of all odd-order coefficients are given by multiplying the signs of the related 
first-order coefficients 

Example Considering the Tables 2 1 and 2 2 again, m the ON-cube 10x0 the variable 
Xi is in affirmation, while the variables X2 and X4 are in negation Therefore, the sign of 
the partial spectral coefficient si is positive and the signs of partial spectral coefficients S2 
and S4 are negative 

The sign of the even-order(2) partial spectral coefficient su of cube 10x0 is positive, 
since the sign is determined by the related first-order coefficients, Si and S2, times -1, 2 e 
(-l)xlx(-l)=l 
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The sign of the odd“Order(3) partial spectral coefficient S 124 of the same cube 10x0 
IS positive, since it is determined by the related first-order coefficients, si ,53 and S 4 t e 
lx(~l)x(-l) = l 

The algorithm is as follows 

Algorithm Calculation of spectral coefficients for completely and incompletely speci- 
fied Boolean functions 

1 For each ON- and DC-cube of degree m, calculate the value and sign of the contribu- 
tion of this cube dto the full n-space spectrum according to the properties described 
previously 

2 The values of all spectral coefficients 6/, except are equal to the sum of all the 
contributions to the spectral coefficients from all ON- and DC- disjoint cubes 

3 For a completely specified Boolean function the value of the dc spectral coefficient 

is equal to the sum of all the partial spectral coefficients corresponding to all of the 
disjoint ON— cubes describing the function, plus the correction factor — 1) x 2”', 
where k is the number of disjoint cubes m the bst of ON-cubes 

4 For an incompletely specified Boolean function the value of the dc spectral coefficient 
So IS equal to the sum of all the partial spectral coefficients corresponding to all of 
the disjoint ON- and DC-cubes describing the function, plus the correction factor 
-(k - 1) X 2^ - ly X 2^“^, where k IS the number of disjoint ON-cubes and w is the 
number of disjoint DC-cubes 

The correction factor — (fc — 1) x 2" compensates for the fact that the cubes over the 
complete n-space have been added k times during the calculation of the k partial spectral 
coefficients A similar explanation applies to DC-cubes as well 

The main advantage of this algorithm is that it can calculate each coefficient separately 
or in parallel So if some of the 2" spectral coefficients are not needed for a particular 
apphcation, then a reduced number of operations can be performed 

Example An example of the calculation of the S spectrum for the four-variable 
incompletely specified Boolean function is shown in Tables 2 1 and 2 2 Fig 2 4(f) 
shows the input data for the algorithm of this section The values and signs of all the 
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partial spectral coefficients for this function has been determined by the properties 1-7 
The results are shown for two cubes 

The spectral coefficients of the first ON~cube in Tables 2 1 and 2 2 (cube 10a:0 of 
degree m = 3) are as follows 

1 within its own m-space, treated as a single minterm; 

So = 6, Si = 2, S2 = -2, S3 = 0 
S4 = ~2, S 12 = 2,si 3 = 0,si4 = 2 
^23 = 0,524 = “2,S34 = 0,Si23 = 0 
^124 “ 2, Si34 = 0, 5234 = 0, 5i234 = 0 

2 Within the full n-space of Boolean function ^(partial spectral coefficients), 

So = 12, Si = 4, S2 = ~4, 54 = ~4 
^12 = 4, 5i4 = 4, S24 = —4, 5124 = 4 

All other coefficients for this cube are zero On the other hand, the spectral coefficients of 
the DC-cube (cube 0000 of degree m = 4, z e , single minterm) are as follows 

1 within its own m-space, treated as a single minterm, 

Sq — 7, 5i — 1,52 — 1? “^3 “ 1 

S 4 = -l,Si 2 = — 1 , 5 i 3 = “Ij’Sh = —1 
^23 = -1,^24 = -1,^34 = '-l,'Sl23 = ”1 
Si24 = —19^134 = -1,^234 = “1,^1234 = “I 

2 Within the full n-space of ^(partial spectral coefficients) the same as within its own 
m-space since it is a single minterm 

In order to obtain the values of all of the spectral coefficients of the whole function, 
except So, the columns of partial spectral coefficients corresponding to all cubes describing 
the function are added (step 2 of the algorithm) The value of sq is obtained by the addition 
of all partial spectral coefficients with the correction factor(step 4) Step 3 is not performed 
because it is an incompletely specified function 

The resulting spectrum is shown in the bottom row of Tables 2 1 ^and 2 2 and, as 
can be easily checked, is exactly the same as obtained by matrix multiplication 
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2.4.1 Spectral CoefRcients for Systems of Boolean Functions 

The algorithms from the previous sections can be modified easily to calculate Walsh spectra 
of systems of Boolean functions In this section, we describe a metl od to calculate spectra 
of a system of incompletely specified Boolean functions 

Let us assume that the functions in the system are in the order ^[1], ^[2], 

9 ^[^]'> where b is the number of functions in the system and the function y^[b] is on the 
rightmost position of the system Then, for a system of incompletely specified functions, 
the following properties hold 

1 The contribution of the spectrum of the function z = 1,2, ,.6 to the total 

spectrum of a system of Boolean functions Stot is equal to the spectrum of the 
function !F[t\ calculated by the algorithm of the previous section, which in turn has 
to be modified by Equation 2 27 

2 The total spectrum of a system of Boolean functions Stot is equal to the sum of all 
the modified spectra of all the Boolean functions m the system 

The contribution of the spectrum of the tih function J^[i] to the total spectrum of a 
system of 6 Boolean functions is denoted in Equation 2 27 by and the spectrum of 
the ith function calculated by the algorithm of the previous section is denoted by Sj[^] 

= 2*’-* X (2 27) 

Conversion of S-spectra to H-spectra for a system of Boolean functions The individ- 
ual first and higher order coefficient values in S and R for a system of Boolean functions are 
related hnearly as in Equation 2 13 However, for the DC-coefRcient(so), the relationship 
is given by, 

= (2 28) 

Example As an example we consider the two functions of section 2 2 2 (page 11) The 
steps of the calculations are shown in Table 2 3 First, the S-spectra of the individual 
functions are calculated using the algorithm of the previous section and are shown m the 
first two rows of the table The modified value of the spectrum of function ^[1] is shown 
m the third row of the table For the function E[2] the modified value of spectrum is equal 
to the original one Then using Equation 2 27 the combined S-spectra calculated and is 
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shown in the fourth row This combined S-spectra is then converted to R-spectra, which 
IS shown in the last row of the table 



So 
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0 

0 

0 

0 

8 

^[1] 

0 

0 

0 

0 

0 

0 

0 

16 

Stot 

0 

4 

4 

4 

0 

0 

0 

12 

Rtot 

12 

-2 

-2 

-2 

0 

0 

0 

-6 


Table 2 3 Table showing calculation steps of spectral coefficients the two functions 
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Chapter 3 


PLA DECOMPOSITION 


Programmable Logic Arrays(PLAs) are important subsystems in the design of digital inte- 
grated circuits A PLA provides a simple and regular layout strategy for Boolean equations 
expressed in two level canonical form, and is usually used to implement random logic (ran- 
dom in the sense that the designer sees no regular structure in the Boolean equations) 
Typical examples are the controllogic for a RISC^ or the control logib« for the microse- 
quencer of a microcoded machine With the addition of latches for feedback, PLAs are 
often used for the combinatorial logic in a finite state machine The optinaization of PLAs 
IS a useful application of computer aided design to the automated synthesis of custom VLSI 
designs 

Techniques for optimizing the structure of a PLA have become weU understood The 
optimization goals are to minimize the area occupied by the PLA and to minimize the 
delay through the PLA The regular structure of a PLA means that the area of the PLA 
IS simply proportional to the number of product terms m the array, and to a first order 
approximation, the delay through the PLA is also proportional to the number of product 
terms (i e , independent of the structure of each product term) 

A PLA macrocell can use input decoders which, group the input signals into pairs [30], 
[31] When the inverters of a standard PLA is replaced with decoders, it is called the PLA 
with decoders The standard PLA is a special case of a PLA with decodersj' being a PLA 
with one-bit decoder The four decodes for each pair are used in the core of the PLA rather 
than the signals and their complements An example of a PLA with input decoders is shown 
in Fig 3 1 A PLA which uses input decoders can always be built with no more rows than 
required by the normal PLA structure [32] However, the number of bteral hues (columns) 
IS equal to that of a standard PLA So, the array size of a PLA with two bit decoders never 
exceeds that of a standard PLA Table 3 1 shows the number of rows needed to realise 
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various classes of functions [32] 
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Figure 3 1 PLA using input decoders 



Standard PLA 

PLA with 
two-bit decoders 

Arbitrary function 
(worst case) 

2n-l 

2n-2 

Syminotnc function 
(worst casej 

2n-l 

3 {n- 2)/2 

Parity function 
worst case) 

2n-l 

2 (n- 2)/2 

Random function of 
lO-variables(average) 

103 

120 


Table 3 1 Number of columns to realize n- variable functions by PLAs (n is even) 

In tins chapter wo will discuss the linear decomposition of PLAs which involves sep- 
aration of the linear (that can be implemented with EX-0 II gates alone) and nonlinear 
(which requires a complete basis such as AND, OR, NOT gates, or NAND gat^s) blocks [7], 
[33] It leads to a cascaded implementation with a linear preprocessor This implementation 
has reduced complexity and is accompanied by an improved testability since the linear part 
IS comjiletely testable for singlc-stuck-at faults and the nonlinear part, iinpleinentcd as a 
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standard PLA has reduced complexity 


3.1 Notations and Definitions 

Let us consider an m-mput, A:-output Boolean function / {0,1}’" {0,1}* The k 

outputs are denoted as /fc-i, /jt- 2 , , /o so that using decimal indices, x'= {xm-i, , x — 

0) € {O,!}”" can be expressed as i = a:t2' and / = ,/o) G {0,1}* by 

f — EJLo The symbols x and / will be used interchangeably to denote the binary 
vectors as well as their decimal representations, for convenience of notation 

The output of / can be described as a vector as in chapter 2, with entries specifying 
the outputs for each input vector a:,0<a:<2™ — 1 

i^=[/(0),/(l), ,/(2’"-l)f, (3 1) 

where f{x) takes on integer values, and 

F^ = [m,m, ,/.(2’"-i)f, (3 2) 

where fi{x), 0 < i < k~l, takes on binary values, and the superscript T indicates transpose 
of the vector 

Definition The autocorrelation of / is defined for its individual outputs /, is given 

by, 



2’"-l 

x=0 

(3 3) 


5 (/.,/.)( 7 .), 0 < 7 - < 2 ’" - 1 

(3 4) 


«=o 


5(/.,/.)(r) can also be computed by the application of the Weiner- Khtnchin theorem [7] 
applying the Walsh-Hadamard transform twice 

= 2-’"lT(17(/.))(r),0 < tau < 2’" - 1 (3 5) 

Where, W(f) is the Walsh-Hadamard transform of / 

The autocorrelation is a measure of the similarity of function outputs at a certain 
distance apart in the Boolean space For example, ^ for all r with Hamming 

weight 1, IS a measure of the number of true minterms at Hamming distance 1, and hence 
can be considered as the closeness of I’s in a Karnaugh map of the function ( Hamming 
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weiglit of a vertex is the niinibor of I’s in the vertex and Hamming distance between two 
vcilices IS tlio inimhor of bits in whidi they diller from each other) 

The definition of autocorrelation can be used to formulate complexity of a function by 
substituting r = 1,2,4, ,2^”^ (of Hamming distance 1), because these are the number 

of true minterms at distance 1 Thus, 

V>(/)=: (3 6) 

t=1,2,4, ,2'"-' 

IS a measure of the function’s simplicity The reason for using this criterion is that two 
minterms that are at unit distance can be combined to form an imphcant, and hence, 
amounts to a reduction in the literal count 

3.2 Linear Decomposition 

The optimal linear decomposition problem is stated as follows 

Find a logic block a, consisting of EX-OR gates, that will minimize the complexity 
measure of the block fa (Fig 3 2) to be implemented over a complete basts 
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Figure 3 2 (a) Oiiginal, (b) Decomposed implementation of a function 

The effect of the EX-OR block is a simple basis translation, a rotation of the Boolean 
space in which / is defined, so as to maximize the criterion ip in the new orientation The 
preprocessor c 7 and the core function fa can be constructed by the following procedure [5] 
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Procedure Let T be an m x m matrix with columns t, so that T = [ro,Ti, . ,r^_i] 
(r, IS m X 1), where r is found as follows 

1 = max(B^^’J\T)), r ^ 0 

2 Let Lq = {0, To} 

Ti is then obtained as follows 

= max{B^^'^\T)),T 6 Lq 

3 Once To,ti have been found, t* can be determined as, 

B(Lf)T^ - max{B^Lf)^^'^'j^j. ^ 

where, 

2S-1 

Ls^l = {^ 5 - 2 , {{Ls^2} © ^s~l)} 

t=0 

Here, © implies bitwise EX-OR of with each element of Ls ^2 

4 a = r""^, where the matrix inversion is carried out mod-2 over the binary numbers 

{ 0 , 1 } 

Once the matrix a is determined, the linearpart can easily be implemented to convert 
any given vector x to by inspecting the rows of a A direct implementation of the 
linear part is given by <%<m Obviously, each EX-OR gate 

implementing a requires as many inputs as there are I’s m the tth row of a In the 
worst case, the matrix a will contain m{m - 1) nonzero entries Hence, the complexity of 
implementing the linear part grows proportional to the square of the number of variables 
m the worst case 

The nonlinear part can be computed as follows Multiply every minterm -,^ 0 ) 

by the matrix a This new vector is a minterm for the nonhnear part 

It is easy to see how this procedure is able to reduce the logic complexity of a function 
The logic simplicity rp was defined earher on the basis of autocorrelation If the maximal 
autocorrelation coefficients correspond to coordinates with Hamming weight 1, then the 
function has relatively low logic complexity because a large number of the minterms are at 
Hamming distance 1 If this is not the case then shifting these coefficients ^to coordinates 
with unit Hamming distance will produce a new function with low complexity which can 
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be generated from the original function by a basis translation The new function (with 
reduced complexity) is f^r and the translation is done by the matrix cr Whether such 
a translation will result in substantial improvements in logic complexity depends on the 
nature of the function itself If all the autocorrelation coefficients are of the same order 
of magnitude than the function outputs do not exhibit any significant correlation Ran- 
domly generated functions and many types of control functions belong to this category 
On the other hand, arithmetic functions, error correcting logic, and many signal processing 
functions have highly correlated outputs and seem to benefit from such a decomposition 
Example Let us consider a two-bit adder fix 3 ,X 2 ,xi,xo) = 2 ( 0:3 -h xi) + X 2 + xq 
According to Equation 3 1 , F = [0, 1,2,3, 1,2, 3, 4, 2, 3, 4, 5, 3, 4, 5, 6 ]^, where the ith row 
of F is the decimal notation og the 3-bit binary output for the adder for the ith row of the 
truth table 

The autocorrelation is computed by Equation 3 4 , 

= [22, 8, 10, 6, 8, 16, 6, 14, 10, 6, 18, 4, 6, 14, 4, 12f 


Using the procedure, after discarding the maximum over the rest of the 

rows IS 18, corresponding to row number tq = 10 (decimal) = 1010 (binary). Lo = { 0 , 10 } 
Without considering the rows of that correspond to members of Lq, the next highest 
correlation coefficient is 16, m row 5 Hence, ri = 5 = 0101, ii = {0, 10,5, 15} 

Similarly, T 2 = 7 or 13 We can choose any one Choosing 13 (1101) arbitrarily, 
L 2 = {0,10,5,15, 13,7,8,2} X 2 remains same whether we choose 13 or 7 because the set 
L 2 contains the linear combinations of both 13 and 7 Similarly, rs = 1 = 0001 Then, 


10 10 
0 110 
10 0 0 
0 111 


To Ti T2 r3 


0 0 10 
1110 
“1010 
0 10 1 

The block implementing a (Fig 3 3) translates x € {0, 1}^ to 2 € {0, 1} The truth table 
for / can be translated with respect to a as follows 
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fcix) = f{Tz) 


For example, let x = (1101) = 13 

'OOlOiri] [o' 

Tz — 

10 10 0 - 1 

.010lJ[lJ [0 

so that /(llOl) = /<,(0010) = ^(2) = 4(= 100) 

iv = [0)1 >4,3,2, 1,2, 3, 4, 5, 4, 3, 6, 5, 2, 3]^ The function is then minimized using 



Non ilntar block as m PLA 


Figure 3 3 Linearized implementation of a 2-bit adder 
ESPRESSO [20] and the cover is given below 
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The implementation of / as a serial combination of a and is given in Fig 3.3 The 
number of literals in the two level SOPE of / is 43 and for fa it is 15 The number of 
product terms for / and fa are 11 and 5, respectively 






The above procedure finds the linear decomposition for an exhaustive list of the out- 
puts of the function However, for an incompletely specified Boolean function the autocor- 
relation can be found out by using Equation 3 5 and the calculation of Walsh-Hadamard 
transform can he done by using the algorithms developed in the previous chapter The 
overhead for the serial implementation is the number of EX- OR gates used to implement o*, 
which may also be implemented by PLAs However, the reduction in literal count, however, 
IS more direct and results in an improved multilevel implementation [34] 

3-3 Results 

The procedure described m the previous section hais been implemented as a computer pro- 
gram “LINDEC”, written in ‘C’ The elfect of serial linearization on some commonly used 
functions is illustrated in Table 3 2 All the PLAs in the table were minimized using 
ESPRESSO In case of the decoded PL A the pairs were chosen by the ESPRESSO ”-do 
pair” option It was found that functions that belong to the class of error correcting and 
translating logic such as code converters, adders and functions with imbedded additions 
benefit most from this approach Symmetric functions which are usually difficult to imple- 
ment efficiently by PLAs also benefit from bneanzation Randomly defined functions are 
not considerably affected by linear decomposition 
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No 

of 

No of rows 

Decoded PLA 

After linear 
decomposition 

Function 

Inputs 

Outputs 

Before 

mimini- 

zation 

After 

minimi- 

zation 

No 

of 

rows 

No of 
decoders 

No 

of 

rows 

No of 

2-bit 

EX-ORs 

Full adder 

3 

2 

8 

7 

4 

1 

4 

2 

2-bit adder 

4 

3 

16 

11 

5 

2 

5 

3 

3-bit adder 

6 

4 

64 

31 

10 

3 

9 

6 

4-bit adder 

8 

5 

256 

75 

17 

4 

15 

7 

3-bit multiplier 

6 

6 

64 

33 

24 

3 

33 

0 

4-bit multiplier 

8 

8 

256 

127 

89 

4 

123 

1 

Binary to BCD 
code converter 

0 99 

7 

8 

100 

45 

33 

3 

33 

6 

BCD to binary 
code converter 

0 99 

8 

7 

100 

36 

30 

3 

36 

0 


Table 3 2 Effect of serial decomposition on some functions 
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Chapter 4 


FAULT DETECTION IN PLAs 


The design of circuits based on PLAs as multiple-output sum-of-products expressions has 
become very popular [35], since the geometric layout properties of PLAs are useful in VLSI 
designs However, the testing of VLSI circuits m general poses difficult problems due to the 
higher density of the circuits and consequent increase in testing complexity For PLAs in 
particular, the methods of random pattern testing have not been proved effective and the 
complete task of test set generation still remains the mam option 

In classical testing a list of fault-free response of the PLA to the test set is required 
For most practical cases, the large storage requirement of such a list makes such test proce- 
dures very expensive Moreover, the computational cost to generate the test set increases 
exponentially with the circuit size Syndrome testing [36] is a very simple yet effective cod- 
ing which counts the number of logical I’s appearing at the output of a circuit The storage 
problem is solved in this method because only the syndrome of the fault-free circuit has to 
be stored Also, the expensive stage of test generation and fault simulation is ehmmated 
When a fault causes the syndrome to change, the circuit is syndrome testable for that fault 
it is said that the test covers the fault Two-level irredundant circmts are syndrome testable 
for their internal lines and can be made syndrome testable for primary input faults through 
additional hardware in the design 

It has been proved [18] that by using a weighted sum of syndromes of all the outputs, 
single stuck-at-faults, single bridging faults and single cross-point faults for all internal and 
input lines in any implementation of PLA are covered 

4.1 Notations & Definitions 

X 2 , , Xji) IS a two-level and-or function of n inputs, represented in the hst of cubes 

form as described m section 2 3 It is represented as, ^ = pi -f P 2 + + Pti (Fig 4 1) 
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Figure 4 1 A PLA 

where xi,T 2 , ,t„ are input lines, [F]j,[J’] 2 , ,[F]jt are the output lines from the k 
output functions, pup^, ,pt are product term lines 

Definition I The syndiome 5(F) of a Boolean function J^{xi,X 2 , ,x„) is defined 

as 5(F) = n (F)/2", where fF(F) is the weight of F and is given by the number of ON- 
vetices in the list of cubes representation of F The syndrome is a functional property 
Thus, various realizations of the same function will have the same syndrome 

ExampU In the function F(zi,X 2 ) = the number of ON-vertices, i e , the no 

of true mintenns is equal to 1, thus th syndrome 5(F) = 1/4 In general, the syndrome 
of an AND gate of n-inputs is 5 = 2~", for an OR gate of n-inputs 5 = 1 — 2~", for an 
EX-OR gate 5=1/2 etc 

Definition 2 Given a set of k output functions [F]i, [Fjj, , [F]^, define the weighted 
syndrome sum IT55 = ^ 2^=1 where wi, W 2 , are set of weights such that 

W 3 = 2’~\ I < s < k 

Definition 3 If the function F is written as F = x,/o+Xt/i (Shannon decomposition), 
the first order Walsh spectral coefhcient r, is given by the number of Ts in /o minus the 
number of Ts in /i, 1 c , r, = lF(/o) - fE(/i) 

Definition J, If F is wiitten as 

F = f.ij /00 + S,Xjfoi + XtXjfio + x.Xj/n, (4 1) 

the second order spec tral eoelliciciit r is given by r,j = lP(/oo)+^^^(/oi)+^F(/io)+fF(/ii) 



Definition 5 For each first and second order coefficient r, and r,_,, R, and J?,j denote 
the corresponding weighted sum of the coeffiaents over all functions [T\ as 

k 

= (4 2) 

5=1 

k 

ily = (4 3) 

5=1 

for the set of weights U), defined above, where [r,], and [r.^], are r, and r.j for output 

4.2 The Classification of Faults in PLAs 

There are three t>pes of single faults to be considered 

Stuck-at- faults One line, or line segment, stuck at either logical 1 or 0 

Bridging faults Two adjacent lines are shorted together assuming the value of the logical 
AND or OR of their own individual signals (determined by the hardware implemen- 
tation) 

Cross-point faults The loss of contact of a connection or its spurious presence 
Each type is considered separately 

The internal single stuck-at-faults m a PLA may be grouped together as shown in 
Table 4 1 with explanatory remarks This is a detailed enumeration of types of faults Later 
these are grouped as to fault effects Each class is exemplified by referring to the layout of 
Fig 4 1 For clarity all vertical lines are labeled with capital letters, and horizontal hues 
with small letters 

Example Considering type 1, which is a s-a-0 fault between literals, is the fault on the 
segment connecting two adjacent literals in the AND array, e g , the horizontal connection 
between the AND connections at xi and X 2 (lme g in Fig 4 1) Now, let us consider type 
7, s-a-0 between products it is the fault on the segment connecting two adjacent product 
lines in the OK array, ( g , the vertical segment between the OR connections of p2 pz 
(line M in Fig 4 1) 

Cross-point faults, t c , erroneous connections between two crossing lines of the arrays, 
can be seen to be of four types 
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Fault type 

1 ) s-a -0 between bterals (c 9 , bne g/O) 

2 ) s-a -1 between literals (e g , line gfl) 

3) s-a -0 sec input (e g , line A/0, line C/0, line 6 / 1 ) 

4) s-a -1 sec input (e g , bne A/1, bne C/1, bne 6 / 0 ) 

5) s-a-0 product bne {eg, bne pi/0, line h/0) 

6 ) s-a -1 product bne (e g , bne pi/ 1 , line h/1) 

7) s-a-0 between products (e g , line M/0) 

8 ) s-a -1 between products (e g , line M/1) 

Table 4 1 Types of faults illustrated with respect to Fig 41 

Growth fault Loss of AND contact between product bne and secondary input line, where 
one should exist (e g , no contact at p 2 and C in Fig 4 1 ) 

Shrinkage fault AND contact between product bne and secondary input bne, where one 
should not be {eg, contact at pi and C in Fig 4 1 ) 

Appearance fault OR contact between product line and output line, where one should 
not bo (e p , contact at p 2 and J 3 in Fig 4 1 ) 

Disappearance fault Loss of OR contact between product bne and output bne, where 
one should exist {c g , no contact at pi and in Fig 4 1) 

A similar enumeration can be applied to the possible bridging faults 

1 Bridging fault between two secondary inputs after inverter (e g , A and C in Fig 

1 1) 

2 Bridging fault between two product bnes {eg ,p, and Pj in Fig 4 1) 

3 Bridging fault between two output lines in OR array {eg , Ts and Tt m Fig 4.1) 

4 Bridging fault between two primary inputs {e g ,x, and Xj in Fig 4 1) 

Tho bridging will bf either AND bridging or OR budging depending on the technology 

4.3 Testability of Stuck-at-faults 

All, change in the syndrome of a function implies that the total number of mrnterms realised 
has' increased or decreased Clearly, an increase or decreas in the syndrome of one output 
causes a corresponding change in the WSS 
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Fault type 

Increase 

Decrease 

1) s-a-0 between hterals 


y/ 

2) s-a-1 between hterals 

V 


3) s-a-0 sec input 


yj 

4) s-a-1 sec input 

V 


5) s-a-0 product line 


yj 

6) s-a-1 product hne 

V 


7) s-a-0 between products 


yj 

8) s-a-1 between products 

V 



Table 4 2 Summary of s-a-f classes 


A PLA implementation is based on sum-of-products-expressions(SOPE) Any given 
stuck-at-fdult affects one or more product terms, but all are affected in an umdirectional 
manner It is demonstrated tbat any internal stuclc“at"fault causes all the affected syndromes 
to intrease or decrease, there being no possibility of cancellation 

Suppose a literal is stuck at constant 1 or is disconnected from the AND hne This 
IS a growth fault [35], and all affected product terms loose a literal leading to an expansion 
in the number of minterms they cover Similarly, if a literal is incorrectly connected to an 
AND line(shrinkage fault), all affected product terms will gain that hteral and there is a 
reduction in the number of covered minterms If a literal is stuck at constant 0, an entire 
product term dis<ippears(disappcarance fault) and WSS of the function will decrease 

h'jauiplf Let us consider a stuck-at-fault of class 1 of Table 4 1 IfbnepmFig 4 1 
IS stuck-at 0, the whole product term pi is ANDed to a result for all input configurations 
Thus, some minterms previously covered by pt loose their coverage causing a decrease in 
the v,eight of allthose functons sharing pt 

Similarly, let a line A be stuck-at-l(class 4 of Table 4 1) Then all produetterms 
connected to A will double in variable xi that is, thay cover minterms both for xi = 1 
j ^ 0 Thus, an increase in the number of minterms covered happens for all functions 


dependent on Xi 

Tablo 4 2 Rivof. i summary of the stuck-at-fault classes as listed previously, with 
lorrespondiiiB iiurcaM./derroasc on the site of the product terms marked beside them 

The evaiiiiiiation of each type of fault demonstrates that a single fault cannot cause a 
literal to lie l.olli removed and added in the sum-of-products expressions, due to the struc- 
ture of the PLA Itself Thus the corresponding fault effect on the WSS is unidirectional 
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Hence, all internal single stuck-at-faults are covered by the WSS However, as will be seen 
in the next section all stuck-at-faults at the primary inputs are not W^^^-testable 

4.3.1 Primary Input Stuck-at-faults 

If an input line is stuck-at-0 or stuck-at-1, both secondary input lines receive the effect of 
the fault, where on the negative line the fault is also inverted This is analogous to two 
separate faults, one of class 3 and one of class 4 (Table 4 1) Prom Table 4 2 it can be seen 
that one fault increases the individual syndromes while the other decreases them Thus, 
there is a possibility of cancellation 

It IS important to determine when and if this cancellation may occur in order to 
maintain WSS testability for the PLA A general result has been proved [18], which is 

given below 


Theorem 1 r./O and r./l are WSS testable tfand only ifR, jL 0, where R, = w,[r, 

Proof Output [P], can be written as 


[•^]j - ^«[/o], + a:.[/i]. 

(4 4) 

Let mo, = < s < k) Now, 


k 

1^5.9 = ^ m,(mo5 -f- mi^) 

5=1 

(4 5) 

-P. = 

5=1 

(4 6) 

i, 

II 

(4 7) 

Now, let us consider x,/0 kP([.?^],)* (m the presence of the fault) is 2mo, 

and WSS* = 


2 ]Ct=i since [/o], is evaluated twice by the test procedure It follows that WSS* ^ 

1T55' if and only if R, ^ 0 Similarly, for x,fl 

Given this result, only each R, need be checked if nonzero, then the corresponding 
input hue is testable by the weighted sum of syndromes, and no cancellation occurs This 
%erifi< ation that each /f, is nonzero is carried out at the design stage and it is not part of 
the actual test At testing time, only the WSS has to be calculated 

Moreover, the spectral coefficients are such that, all of them are either odd or even 
for any given completely specified function This is because of the orthogonality of 
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the transformation matrix T [2] For the set of weights w„ it is clear that = [r,]j + 2 A 
for some A, since the weights are multiples of 2 (except tui) Thus R-tv even or odd 
according as [r,]i is even or odd, this being determined by whether is of even or odd 
weight However the ordering of output functions can be permutated to ensure that a 
function of odd weight becomes the right most output, i e , This selection is done at 
the design time 

for an incompletely specified function, in the PLA implementation the don’t care 
product terms are not implemented Thus it becomes equivalent to a completely specified 
function, with its DC-cubcs converted to OFF-cubes and the above considerations are 
valid 

I he onl\ tune when the above permutation method faals is when all functions being 
rcah/ed are of even weight After calculating R^ in such a case, and finding it equal to 0, the 
simplest hardware change done is the addition of an extra output line This extra function 
realized in the PLA should have only one minterm, yiy 2 € {x„f, }, 1 < t < n This 
extra output is made such that its wi is 1 In this way the value of R^ becomes +1 or 
— 1, and the PLA becomes IF .9 S' testable for primary input stuck-at-faults 

4.4 Testability of Cross-point Faults 

The four t\pc‘s of cross-point faults arc covered in a similar manner to the stuck-at-faults 

A shrinkage fault betw'c'en a product line p, and a secondary input line causes the 
coverage* of p, to either be halved, if it is not already connected to or to go to zero, 

if it IS already connected to x, Thus the weights of all functions involving Ps are either 
dec reased The same applies for a contact between p, and f, 

Conversely, for a growth-fault causes the Ps involved to become independent of the 
variable z, Thus, the product term doubles in size along x, and the weights for all the 
functions sharing it increa.se 

For an appearance- fault, an extra product term is realised for the function involved 
causing an incrc'asc* iii the number of minterms covered (note that there will be no change 
if thc'rc* IS rc'dunciaiic y ) A disappearance-fault for a function reduces the number of 
minterms covered 

I'hiis, all internal single cross-point faults are covered hy WSS 
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4.5 Testability of Bridging Faults 


The five classes for AND bridging faults are considered here m detail Similar results foUow 
for OR bridging faults 

1 Bridge between x. and z, Any output [J^l that depends on z. can be written as, 

[J"], = Ax, + Bx, + C, (4 8) 

for some A, B,C Then under the AND bridging fault between z. and z„ it becomes 


+ £(a:,z,) + C 

= C (4 9) 

f learK U is less than the fault-free case and 14^5" 5' changes accordingly 
2 Biidqt Ih hiHtn J, and z^ This fault causes both z. and Zj to be replaced by (z,Zj) 

[^], = Az.Zj -t- Rz.Zj + Cz.Zj 4- Dz.Zj + E (4 10) 

for some 4, H, (\ 1), E changes to, 


[^l' = + C(z,Zj) + D, (4 11) 

With a <()rr(‘sp()iKling docreaso m \\'{[calF\s) s-nd WSS 

3 Hi Ih twi i n pimliK t tt nn hru.s and pj This fault causes both the product terms 
to be r(‘pla(<Hi by [pxVj) output = Pi + P2 + +Pt+ +Pj + becomes 

[^]* = Pi + P2 + ■ + ?.?;+ , (4 12) 


With d r(‘sultint^ decrease in weight 

4 Hr id(j( Im tuH ( n output hius [F]^ and \T\ This type of fault m the OR array causes 

tlie two outputs to become 

\n: = irr, = [nm, (■* w) 

With d (ifTrease in weiglits of both the functions and 

5 Hr idgt Ih tw( ( n pi irnaiy inputs^ Xx end Xj This is a special case and we consider it in 

the next sec tion 
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4.5 1 Primary Input Bridging Faults 


In the case of an AM) bridging fault between two primary inputs, x,,Xj, the effect is earned 
to both secondary input lines with the following multiple substitutions taking place 

• each of x, and Xj is replaced by x^Xj 

• each of x, and Xj is replaced by f.fj 

If [J], IS written as liquation 4 10, then under the fault it becomes, 

[/■], = Ax,Xj + Dx, 4- Dxj + F, (4 14) 

and a ( aiiccllat ion potential exists, since there is no guarantee that the number of minterms 
uiufonnh decreases or increases 

'I’o isolate this possIblllt^, a condition similar to the one previously discussed for stuck- 
at faults on pnniars inputs can be prosed If T is expressed as in Equation 4 1, for an 
AM) bridging f.iult hetv,ec-n tv,o priniar> inputs, the subfunctions /oi and fio assume the 
same value <is fim '1 bus no c aiucdlation can occur if and only if 


H’(/oi) + in/io)y^2VF(/oo) (415) 

'I h s condition can he expre-Mui in the spectral domain using the first and second order 
spectral coefficients, and hc-rc' the general result is proved [18] 

Thc'oreni 2 An \\I) biulqing fault (Htw((n two primary inputs and Xj isWSS 

(< '•tahlt if and anhj if U, Hj } 11,^ f 0 


I’rooJ We h.ise. 


[/■], - j,fj[/(X)], 4 x.Xj[/oi], + x,Xj[/io], + x.Xj[/n]s (4 16) 

la^t mo, -- U ([/no],), Till, = »'{[/oiL ), tiij, = H^([/ioU, ^3. = = 

minibei of outputs Now, 


k 

\\ SS -- ^^'•'(”'0’ + + "^25 + '^3s) 

ir 1 

K 

h\ -- 

ir-I 

k 

= ^ U>,(mos + T«l5 - TTl2s - TTlSi) 

,= 1 


(4 17) 
(4 18) 
(4 19) 
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k 



(4 20) 



k 


= - mi, + m 2 , - m 3 ,) 

5=1 

(4 21) 

k 

5=1 

(4 22) 

k 

= tWiCmo, - mi, - m 2 , + m 3 ,) 

(4 23) 


a=l 


With AND bridging between Xt and Xj^ the weights become 

W([foi],) = = M'([/ooU = m„. 

Thus, uja(3mo, + mz^) 

It follows that WSS* ^ WSS if and only if i2, + Rj ^ 0 Similarly, for an OR 
bridging fault, it can be proved [18] that bridging fault between primary inputs x* and Xj 
IS testable if i2, + Rj -Mtj ^ 0 

Thus at the design stage the i2i, Rj and R^j for each pair of physically adjacent input 
lines (there are n — 1 such lines), should be checked If the calculated sum is 0, for any t 
and then the same change in hardware as suggested in the case of stuck- at-faults, % e , 
addition of an extra output line works effectively 

In the case of AND bridging faults note that the extra minterm, previously written as 
2 / 1^2 Vn^ where each is a literal being either x, or x,, is best implemented as XiX 2 Xn 
This particular product term will affect only the total weight of the foo decomposed portions 
of the output functions by adding one to it and leaving the other two total weights, i e , 
kF(/oi) and W{fio) unaltered This remains true for any pair of (xi,Xj) since, it is always 
just W{foQ) that is affected 

This technique is described in the algorithm below In the first step the PLA is 
expressed as a hst of disjoint ON-cubes for each output with the algorithm of section 2 3 
Then, the calculation of WSS becomes easy, since weight of each function is the total 
number of vertices m all the cubes m the list The number of vertices in a cube C is 
given by, W{C) = 2°^, where d = number of don’t cares in C Also the same list is used to 
calculate the 1st and 2nd order spectral coefficients using the algorithm described in section 
2 4 Then it recursively checks for the s-a-f and bridging fault testability at the primary 
inputs If, either for a bridging fault or for both a s-a-f and a bridging fault at the primary 
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inputs, the PLA is found untestable, then an extra output line is added, which implements 
only the 0-th minterm However, if it is untestable for only a s-a-f at primary inputs, the 
output lines are reordered so that corresponds to an output of odd weight making the 
corresponding s, 0 

Algonthm Checks the WSS testability of a PLA. If found untestable, modifies it to 
make it WSS-testable 

for each output, 

express the PLA as a hst of disjoint ON-cubes, 
bdg = saf = 0, 
do { 

for each output i = 1 until k 

calculate the 1st and 2nd order spectral coefficients 
for each i = 1 until n, /* n = no of inputs */ 
if {Rt == 0)saf = 1, 

for each t = 1 until n — 1 /* Check for each adjacent pairs of 

primary inputs */ 

if {{Rt + -Ri+i + 2jR,i2,+i) == 0)bdg = 1, 
if {bdg == 1) 

if product term (a:, = 0, 1 < i < n) is present in the PLA 

add an extra output [jP]j having only this product term, 

else 

add the 0-th minterm as a product term in the PLA, 
add an extra output having only the 0-th minterm, 

else if {saf == 1) 

selectflag — 0, 
do { 

select an output line r randomly, 
if(TP'(Fr) is odd) selectflag = 1, 

}while {select flag ^ 1), 
interchange the outputs [F],. and [F]^ 

} while {saf == 1 or bdg == 1), 
calculate the of the PLA, 

4.6 Coverage of Multiple faults 

The assumption of single faults of any type may not be sufficient for full testability However, 
it may be reasonable to expect that, a set of faults occurmg simultaneously are not totally 
random in their appearance and distribution The testability criteria as discussed in the 
earlier sections of this chapter for single faults by W55 are based on the unidirectionality 
of the change m the syndromes Thus it follows that any combination of faults having the 
same unidirectional effect is fully covered hy WSS 

The result acquires more sigmficance by noticing that the most common faults come 
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from the error m programming the PLA, resulting in too many or too few conn*ections, or in 
shorting of a set of adjacent lines Hence the WSS is a strong testability condition for the 
largest and most common set of multiple faults, which indeed lead to unidirectional effects 
[18] 

4.7 Implementation and Results 

This algorithm has been implemented as a computer program, ‘‘SPECTEST”, written in 
^C’ It has the option of calculating the WSS only, checking for complete kPiy^-testability 
and making the necessary hardware modification, and finding the fault in a faulty PLA if 
the fault-free PLA or its syndrome is given as input It has been tested on a set of MCNC 
PLA benchmarks, and some commonly used functions The results are shown in Table 4 3 


Example 

Testability 


Testing 

time(s) 

S-a-f 

Bridging 

5xpl 

V 

V 

55961 

0 56 

Z5xpl 

V 


40768 

142 

9sym 

X 


849 * 

2 17 

Z9syin 

X 

^/ 

841 * 

4 63 

rd53 

V 

V 

76 

0 24 

rd73 

V' 

\/ 

449 

0 98 

rd84 

V 

v/ 

1636 

2 79 

sao2 

V 

^/ 

1417 

1 03 

add2 

V 

V 

48 

0 11 

add3 

V 

V 

448 

0 34 

add4 

V 

y 

3840 

2 98 

mult3 

y 

y/ 

784 

0 28 

mult4 

V 


14400 

2 93 


Table 4 3 Table of test results of PLA benchmarks implies testability and x implies 
untestability ★ implies it is the WSS of the modified PLA The times are given for s-a-0 
faults at primary inputs 
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Chapter 5 


Conclusion and Future Work 


The essential relationships between classical and spectral methods used in the design of 
digital circuits have been stated Based on these relations, algorithms for generation of 
spectral coefficients for both S and R spectra for completely as well as incompletely specified 
Boolean functions have been deshown Such detailed interpretations of S and R spectra 
are important not only from the point of view analysis and synthesis of digital systems, but 
also for the generation of test sets and design for testability 

Linear decomposition, described in chapter 3, has no equivalent Boolean domain 
solutions [5] It is a good supplement to other synthesis techniques because irreduces the 
logic complexity of many functions that are usually difficult to synthesize efficiently in the 
Boolean domain It has also been found to improve testability of implementations because 
EX-OR networks are known to be easily testable for single stuck-at-faults Linearized 
implementations are also known to be syndrome testable, which makes them suitable for 
built-in test 

The weighted syndrome sum for all outputs in a PLA described in chapter 4 is a 
sufficient test to cover single stuck-at-faults, bridging faults and cross point faults The 
model adopted did not exphcitly consider all multiple faults However, all multiple faults 
belonging to classes possessing the same direction of change in their WSS are fully cohered 
The weighted sum of syndromes is a very simple yet very powerful coding for the testing 
of PLAs and its straightforward implementation makes it suitable both for easy design for 

testability (DPT) and built in self test (BIST) 

Future work in this area can be development of new decomposition methods for 
systems of incompletely specified Boolean functions based on the representation of the 
Hadamard- Walsh spectrum presented The properties of such decompositions make them 
very suitable for design using Field Programmable Gate Arrays (FPGAs) Also, some work 
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on the area of multiple fault detection using the syndrome technique can be taken 
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